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1 Introduction 

The multidisciplinary field of spintronics has attracted keen scientific interest 
from the viewpoints of both fundamental and applied physics. [T] The key 
notion there is how to carry the intrinsic spin magnetization in a well defined 
fashion. [21 El H] The current induced motion of a magnetic domain involves 
a torque transfer process from the conduction electrons which generates a 
sliding mode of the domain wall, i.e. a motion without distortion and without 
tilting and twisting [51 [6]. Naturally speaking, the spin current is understood 
as the deviation of the spin projection from its equilibrium value. 

Quite recently, we proposed a new way to generate a spin current in 
the so-called chiral helimagnets where the axial coupling between the neigh- 
boring spins called Dzyaloshinskii-Morya (DM) coupling fixes the left- or 
right-handedness of the helical spin arrangement over macroscopic scales. [7] 
By applying the static magnetic field perpendicular to the helical axis, the 
magnetic kink crystal (chiral soliton lattice) is formed [H [9] and the inertial 
motion of the kinks (Galilean boost) triggers of the transport of magnon 
density (spin current) inside the magnetic background [7\. Starting with the 
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effective lattice Hamiltonian and its continuum version, we wrote down the 
Lagrangian including the appropriate Berry phase term and introduced the 
sliding solution of the soliton lattice by inspection. Then, we constructed the 
canonical Hamiltonian for the collective coordinate variables which describe 
the center-of-mass motion of the whole soliton lattice, based on the Dirac's 
prescription for the singular Lagrangian. After these procedures, we suc- 
ceeded in describing the internal motion of the kink crystal state. The most 
important notion is that the "spin phase" directly comes up in the observable 
effects through the soliton lattice formation. In our mechanism, the current 
is carried by the moving magnetic kink crystal, where the linear momentum 
has a form, P = 2ttSQ + MX. The topological magnetic charge, SQ, merely 
enters the equilibrium background momentum 2irSQ, while the collective 
translation of the kinks with the velocity X gives the mass M. 

From a mathematical viewpoint, however, the starting Lagrangian with 
the Berry phase and the axial coupling [see Eq. ([3]) below] is not manifestly 
Galilean invariant. Then, the following problems remain highly non-trivial: 
(1) how to describe the stability of the ground state under the Galilean boost, 
and (2) how to justify the existence of the linear momentum as a conserved 
Noether current. To guarantee the stability, Lie point group symmetries 
of nonlinear differential equations (DE) describing dynamics of the system 
should embrace the whole transformation. To justify the existence of the 
conserved current, we need to establish appropriate conservation laws from 
the analyses of variational symmetries. The latter plays an important role in 
the study of magnon transport. In this paper, to clarify the above mentioned 
problems (1) and (2), we start with the continuum version of the chiral XY 
model under the magnetic field that describes the chiral helimagnet with 
soliton lattice. Then, we apply the Lie group methods and symmetry analysis 
to the model. Although the model and the method of the present analysis 
have been well established, we are definitely motivated by how to justify 
the emergence of the spin current in real chiral helimagnets that serve as 
potentially promising devices in spintronics. 

The paper is organized as follows. We discuss the model of chiral heli- 
magnet in Sec. I. In Sec. II we consider the problem of finding the Lie point 
symmetries of the DE system of the continuum theory. The main question 
studied here is how to find group-invariant solutions of the differential equa- 
tions. The variational symmetries are deduced in Sec. II. and we establish 
here explicit formulae for the conserved densities and density currents in- 
volved in the conservation laws. Conclusions are given in Sec. IV. 
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2 Model 



To describe chiral spin ordering that occurs along a crystallographic axis, say 
z-axis, the relevant t model Hamiltonian is given by 



H 



-J SA + i -V-J2 [3 x - 2fi h (1) 



The ferromagnetic interaction with the exchange coupling J (first term), and 
the parity-violating DM interaction, described by the uniform axial vector V, 
between the nearest-neighbor spins lead to the long-period incommensurate 
helimagnetic structure with either left-handed or right-handed chirality. The 
direction of the mono-axial T> vector is fixed by the crystallographic symme- 
try that we do not get into the detail[7j. The third term is the Zeeman term 
with a transverse magnetic field with its strength being h. 

In a continuum approximation, we obtain the stationary conditions for 
the 9 and ip variables that lead to the coupled nonlinear partial differential 
equations given by 

S 6 t + JS 2 (sin 9 <p zz + 2 cos 9 9 z tp z ) - 2VS 2 cos 9 9 Z - 2fx hS sin (p = 0, 

S sin 9 ip t +J'S 2 (sin 9 cos 9ip\ — 9 ZZ ) —2VS 2 sin 9 cos 9cp z —2fj, hS cos 9 cos if = 0. 

(2) 

The energy functional is then given by 

H = J dzly [9\ + sin 2 9ip 2 z ] - VS 2 sin 2 9cp z - 2^ hS sin 9 cos if j . 

To describe the dynamics in an appropriate way, we need the Lagrangian 
density, 

L = S(cos9-l) Vt -^[9 2 z + S m 2 9 V 2 z ] 

+VS 2 sin 2 9ip z + 2fi hS sm9 cosip (3) 

where the first term is the Berry phase term. The ground state configuration 
is then the kink crystal (chiral soliton lattice) state specified by the stationary 
solution, 

V9(z) = 2cos~ 1 sn{—z,Kj , = 0, (4) 
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where the helical pitch in the zero field is given by q = T>/J, m = ^2/i h/ 'JS, 
and sn is the Jacobi elliptic function with the elliptic modulus k (0 < 
k < 1) determined through minimization of energy per unit length, n/m = 
4:E(n)/irqo, with E(k) being the elliptic integrals of the second kind.|8] In. 
Ref. [7], we constructed the sliding solution of the kink crystal by intro- 
ducing the position of the kink center Z as a dynamical variable. Based 
on the collective coordinate method and the Dirac's canonical formulation 
for the singular Lagrangian system, we derived the closed formulae for the 
mass, spin current and induced magnetic dipole moment accompanied with 
the kink crystal motion. 

3 Classical Lie symmetries and Galilean boosts 

Now, we perform the classical Lie symmetry analysis to manifest the existence 
of the Galilei-boosted solution of the stationary solution (J4]). The system of 
the second order DEs, ([2]), are eppressed as F 1 = 0, F 2 = 0, where 

F 1 = S sin 9 9 t +JS 2 sin 6 (sin 9ip zz + 2 cos 9 9 z cp z )-VS 2 sin 29 9 z -2fi Q hS sin 9 sin <p, 

F 2 = S sin 9 ip t +JS 2 (sin 6* cos 6^5; — 9 ZZ ) — 2VS 2 sin 9 cos 9ip z — 2 fi hS cos 9 cosip, 

(5) 

with two dependent variables u = (9,(p), and two independent variables 
w = (z,t). The first equation in (jSJ) is multiplied by the factor sin 9 that 
is necessary for the further symmetry analysis. The symmetry condition 
requires that the system must hold in the transformed variables whenever it 
holds in the original variables. 

In order to take into account the derivative terms involved in Eq.fjSJ), the 
infinitesimal generator is the prolongation introduced as [13] (i = 1 means 
"2", i = 2 does "t") 




(*<j) (i<i) 



where the first prolongation is given by 

C, = A (C) - g Di (0 - 9 t Di (77) , Xl = A (X) - (PgDi (0 - <p t Di (v),(i = *, *) 
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and the second one is determined from the first ones 

Cy = D 3 (d) - e zz D 3 (0 - 0^, (17) , = ^ (Xi) - ^ (0 " ^ fo) , 

C 2j = I>i (Ca) - zt D 3 (0 - 0^ (17) , X 2J = Dj (Xa) - ^ (0 - ^ fo) • 

Here, the Z?j are the total derivatives with respect to z and t relevant for the 
system (J5J). 

The system of DEs is invariant under the group G if 

XF i | J r = 0, i = l,2. (6) 

The invariance condition provides a partial differential equations involving 
the unknown infinitesimal functions with their derivatives and products of 
the partial derivatives of the dependent variables. After replacing 6 ZZ and ip zz 
from Eq.([5]) and splitting the resulting equation into powers of derivatives of 
the dependent variables (shown below at left), handled as independent ones, 



the overdetermined system XFi\ F = is obtained 

sin 2 07^ = (7) 

<p zt z :sw 2 0r lg = (8) 

<p gt <p z :am 2 6ri lp = (9) 

9 t 9 t :£ e an9 = (10) 

9 m :^sin9 = (11) 

9 t : -( cos 9 - r/ t sin 9 + 2£ z sin 9 - x^ sin 9 + ( e sin = (12) 

<p t :( v sm9 + x e sm 3 6 = (13) 

0*<Pz ■ ~Xe v + 2C + Xe v cos 20 - C e sin 26 = (14) 

2 :2 Xe cos0 + x ee sin0 = O (15) 

ifl : x w sin 2 6 + x e cos sin 3 + sin 20 = (16) 
9 Z : - JSxe, cos 20 - f t sin + JSXfe + 2£>SC + JSx, sin 20 

-P5sin20(^-^ + C e ) =0 (17) 



cp z : -JS sin 2 9£ ZZ +2JS sin 2 0x^-2£>S cos sin 3 9x e +JS sin 29( Z -VS sin 20C V = 

(18) 
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1 : ( t sin 9 — 2fi h cos (p sin Q\ + JS sin 2 9x zz — 2fx h sin 2 9 cos 9 cos 

— "DS 1 sin 29( z + 2/x /i cos 9 sin — 4// /i sin 9 sin </?£ 2 + 2/i h sin 6> sin <£>x v = 0. 

(19) 

The overdetermined system XF 2 \p = is written as 

2 

9 zt :V Z = (20) 

9 zt 9 z : Ve = (21) 

0*t¥>* : ^ = (22) 

ip z ip t :^sm9 = (23) 

^0 t :&sin0 = O (24) 

<Pt : C cos 9 - rj t sin 9 + 2£ z sin 9 + x^ sin 9 - ( 9 sin 6» = (25) 

^:g + X,sin 2 # = (26) 

0z<P z ■ Ce v - C v cot 9-Xe cos 9 sin = (27) 

^ : Cw = (28) 

V? 2 : 2C W - 2C cos 2£ + (£ g - 2 Xv ) sin 2£ = (29) 
ip z : -2JSC VZ - 2VS cos 2^C - £ t sin - VS sin 20£ 2 + sin 29 Xz 

-VSsm29x v + VSsm29( = O (30) 

2 : J (£ zz - 2C e J - 2V cot - V sin 20 Xfl = (31) 
1 : — JS( ZZ — Afi h cos 9 cos (p^ z +2 fx h cos 9 cos (pC d +Xt sm9+2fi Q hsin9 cos(f( 

- VS sin 29 Xz + 2^h cos 9 sin VX - 2 f i h^( ip = (32) 

Sill C7 

The general solution of the determining equations is presented in Ap- 
pendix A. The final result is 

Z = C U r) = C 2 , x = 0, C = 0, (33) 

where C\^ 2 are the constants. Therefore, every infinitesimal generator of Lie 
point symmetries is a linear combination of X 1 = d z and X 2 = d t 

X = dd z + C 2 d t . (34) 
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We seek invariant solutions of the Eqs.flS]) under symmetries generated by 
The characteristic equations are 



dz dt d8 dip 

c[ = c~ 2 = T = T' 

which have three functionally independent first integrals 



2i = z — vt, 1 2 — 9,Z 3 — y?, 



where v = C\jCi is a velocity, and Xli — (i — 1, 2 
We can insert first integrals having the form 



3). 



9 — 9 [z — vt) , ip = (p[z — vt] 



(35) 



into the original system (j3J). This yields the simplified partial DE system 
after introducing the new coordinate z — z — vt. Thus, the invariant solution 
includes the sliding mode ( |35|) which can be found by inspection, but we clas- 
sify all invariant solutions, using the structure of the Lie algebra. Therefore, 
we proved that the Galilei-boosted solution certainly exists. 

4 Variational symmetries and conservation laws 

Next, we construct the conserved Noether current to justify the existence 
of the linear momentum that eventually carries the spin current in chiral 
helimagnets. A derivation of conservation laws is based on so-called Noether 
identity [Hj 



Y (L) + A • (?L) = W a E a (L) + DiN*{L), (a = 1, 2; % = z, t) (36) 



The dot in Eq. (|36|) means the differentiation rule Di ■ \CL) = Di (C) L + 



where the canonical vector field eligible for the Lagrangian ([3]) 



Y = W 1 d 9 + W% + D z (W 1 ) d 9z + D z (W 2 ) d Vz + D t (W 2 ) d, 



with the characteristics 



w 1 = v 1 - ee t - eo z , w 2 = v 2 - eVt - ev,- 
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In Eq. ( 1361) N l are the Noether operators given by the expressions 

d 



N l = C+ (V a - e<) ( + £ (-l) s D n . . . D 3t 

\ 1 S>1 

+ j2 (-i)- Dkl ... Dkr if - eut) + £ 

r >l \ %-'r S >1 



UU ijl-js 



■l) s D jl ...D j 



ua ik 1 ...k r j 1 ...j s 
(37) 

and Euler-Lagrange operators are 

The variational symmetries on the Euler-Lagrange equations may be found 
via the Noether identity 

Y(L) + D s .(eL)+D t .(?L)=0. 

In explicit form it is written as 

o = n l L +n 2 L v +{p z (w 1 ) + eo zt + ?o zs ) l 6z +(d z (w 2 ) + £ Vzt + £>J L Vt 

+ (a (w 2 ) + ev« + evj ^ + (a (D + a (^)) l 

with the derivatives 

JS 2 

Lg = —S sin 9(p t — sin 29ip 2 z + VS 2 sin 29 <p z + 2fi hS cos cos <p, 

= — 2 / u /iS'sin6 l sin<y2, L 8z = —JS 2 9 Z , 
L Vz = -JS 2 sin 2 Q^l + VS 2 sin 2 0, L Vt = S (cos - 1) . 
After a little manipulation one gets the following equation for the unknowns 

= r/ 1 I —5 sin 9 (p t — sin 29ip 2 z + VS 2 sin 29 tp z + 2[i hS cos 9 cos (p j 

—T] 2 2fi Q hS sin 9 sin </? 
-JS 2 2 $ + 9 zV l e + tptf - 9 t [C + 9£ + Vz Q - 9 Z [C + 9 z C e + 
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- ( JS 2 sin 2 6<p z - VS 2 sin 2 6) 

x M + ^ + - n K + ^ + <o4J - ^ Si + + ^1 ) 
+s (cose - 1) (t? 2 + otfo + ^ - <p t [$ + e£ e + vtQ - ^ z [it + Ot& + vtQ) 



The equation is solved as an algebraic equation with respect to the partial 
derivatives of the dependent variables, handled as independent variables. 
The equation leads to a determining equations with respect to the unknown 
infinitesimal functions £ 2 , £ , 77 1 and if. The details of the calculations are 
relegated into Appendix B. As a result, we obtain, C% t 2 are constants, 

? = C U C = C 2 , V l = V 2 = 0. (38) 

Using the found variational symmetries the operators N l are firstly con- 
structed from formulae (1371) and the corresponding conservation laws are 
then calculated from 

Di [N l L)= 0. 



The Noether operators modify the Lagrangian density into a conserved quan- 
tity C % = N l (L) with a zero total divergence. 

(i) The choice C\ — 1 and C 2 = corresponds to the symmetry group 
under translations in time t — > t + a. Then 



-do t rt ^ t ' ^ ^a^ - 

By introducing the functions 

C* = #*(L) = sin 2 y? 2 - + VS 2 sin 2 5 p z + 2fi hS sin cos <p, 

and 

C 2 = iV 2 (L) = JS 2 9 z e t + JS 2 sin 2 9 <p z (p t - VS 2 sin 2 6 <p t . 
one gets the energy conservation law D t (C) + D z (C z ) = 0, or 

D t (S) = D z (Jf), 
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where 8 = — C l and Jz = C z are the energy density and the density of 
energy current, respectively. 

(ii) The choice C\ — and Ci — 1 yields the variational symmetry under 
translations in space z — > z + a. Then 

n* = -e z — - <p—, n z = i- e z — - ip t 



•do t Tz d^ " z de z ^ z d^ z 

By the same manner one can obtain the momentum conservation law 

A (Pz) = A (T^) , 

with the momentum P z = —C l = —N l (L) 

P z = S(l-cos9)<p z , 
and the canonical energy-momentum tensor T zz = C z = N Z (L) 

JS 2 

T zz = S (cos 6 — 1) cp t H ^— + sin 2 9 Lp 2 ^j + 2fj, hS sin # cos 

Therefore we proved the existence of the linear momentum as a conserved 
Noether current. Then, we are ready to define the longitudinal spin current 
carried by the kink crystal by writing the linear momentum per unit area as 

P Z = S (l-cos9)ip z dz = 2nSQ + Mz, (39) 



where L is the system size. The topological magnetic charge, SQ, merely 
produces the equilibrium background momentum 2irSQ, while the collective 
translation of the kinks with the velocity z gives the inertial mass M of the 
kink crystal [7]. 



5 Conclusions 

In this paper, motivated by the spin current problem in chiral helimagnet, 
we rigorously proved the hidden Galilean invariance embedded in the chiral 
XY model under the magnetic field. The Lie group analysis is applied to 
the differential equations of the continuum theory of the chiral helimagnet 
with the parity-violating Dzyaloshinskii-Morya coupling under a transversal 



10 



magnetic field. Lie point symmetries and the invariant solutions under these 
symmetries are found. They present sliding solutions that come up as a 
consequence of both a breaking of spin rotational symmetry by the external 
magnetic field and a parity violation due to DM interaction. We found that 
variational symmetries are related with translations in space and time, the 
corresponding energy and momentum conservation laws are derived. We 
therefore succeeded in justifying the existence of the transport spin current 
in chiral helimagnet. 

Appendix A 

From Eqs.flTJCEI]) and fl2UH24p one obtain r] = -q(t) and f = £(z,t). From 
Eq.(H5D we get x = Xo(z,t,ip) - Xi(*> *, <p) cot 9 and ( = Ci(z, t,<p)9 + 
£ Q (z,t,<p) from Eq.dSSJ. 

(i) Let us consider x = Xo( z i t,ip). A substitution of the results for x an d 
C into Eq.flT2D yields 

— (d9 + C ) cos 9 — r] t sin 9 + 2£ 2 sin 9 — Xo v sin 9 + d sin 9 = 

and after splitting over the 9 variable we obtain ( = (( = ( x = 0) and 

-Vt + 2£ 2 - Xo v = 0. (40) 

From Eq. (fT3l) we get Xe — that agrees with the choice x — Xo( z it,(f) an d 
transforms (1271) into identity. A sum of Eqs. ffT2l and fl25l) gives % = 2£ 2 , 
and, therefore, £ 22 = and Xo v — (from Eq.[4"0]). i.e. x — Xo( z ^)- Then 
Eqs.( ll6f2l?j) become identities. Eqs. (117l30j) coincide with each other and may 
be written as 

-C t sin 9-VS sin 20£ g + JS sin 29 Xz = 0. 
The relation splits into 

& = o, -x>e, + Jx z = o, 

that gives x zz = 0- Eqs. (pi5B are identically fulfilled. The last Eqs. (11911521 
take the form 

cos ip sin 9x + 2 sin 9 sin 2 = 0, 
— 4/x /i cos 6* cos <^£ 2 + Xt srn $ ~~ sin 20%^ + 2/i /i cos 9 sin = 0, 
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respectively. Since x does not depend on 9 and ip, after splitting over these 
variables we get 

Xt = °> Xz = o, x = o, = o. 

As a result we have eventually £ =const, 77 =const, x = 0, C = 0- 

(ii) Now we take x — ~ Xi( z i c °t 0- We prove that only Xi — 
satisfies the determining equations. Indeed, a substitution of x m to Eq. (l2"51) 
and a splitting over variable gives 

Ci = o, Co = Xhp, Vt = 2 L- 

This means that ( = ( (z,t,<p) and £ 2Z = 0, and Eq. (fl2^) becomes identity. 
From Eq. fT26l) we obtain £ = —\i that results together with ( = Xi v in 

X = — [a(z, t) cos ip + b(z, t) sin <p\ cot 9, 

( = —a(z, t) sin <p + b(z, t) cos (p. 

Eqs.f ll3|fT4|27lll6|29|18ll3ip are identically fulfilled. After splitting over 9 and 
ip variables Eq. (!30|) turns into 

sin#: £ t = 0, sin2#: £ 2 = 0, 

cos ip : Ja z + Vb = 0, sin </? : Jb z — Va = 0, 

cos 2 6* cos <£> : Vb + Ja z = 0, cos 2 # sin ip : —Va + J'b z = 0, 

that yields £ =const, r\ =const, and a z = — (V / J) b, b z = (V / J) a. Eq. flTTI) 
gives the same result. From the last Eq. fll9l) we obtain 

sin^sinyj : a t = 0, sin # cosy? : b t = 0, 

i.e. a = a(z) and 6 = b(z) and 

sin cos 9 cos y> : J7"a 22 + 2T>b z = 0, 

sin 9 cos sin y? : Jb zz — 2T>a z = 0. 

Together with the previous results these equations give a z = and b z = 0, 
therefore a = b = 0, and we get eventually x — C = 0- Then Eq. (!32|) will be 
an identity. Eventually, the result (133|) is again obtained. 
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Appendix B 

A set of equations 



el ■ Co = o 
fa,- £} = o 
o% : el = 

0*0tV> z ■ C = o 

yields immediately £ 2 = £ z (z,t) and £ = £\z,t). A remaining part acquires 
the form 

Mt: £ = (41) 

^: ^-|^ + Kt = (42) 

0^: 77^ + sin 2 0^ = (43) 

<^ 2 : -±8^207/ + sin 2 0^ -sin 2 0r/ 2 - i sin 2 + =° ( 44 ) 

t : (cos0- l)^ 2 , = =>rfi = (45) 

^: -J^ + £>sin 2 0r? 2 = (46) 

(pt : _ sin 0^ + (cos 0-1) (r/ 2 + £) =0 (47) 
^ : PS 2 sin 20 r? 1 - JS 2 sin 2 0r? 2 + VS 2 sin 2 0r/ 2 

-X?S 2 sin 2 6^1 - S (cos - 1) g + VS 2 sin 2 (f « + g) = (48) 



1 : 2fi hS cos cos ipr] 1 — 2fi Q hS sin sin r? 2 + PS* 2 sin 2 r/ 2 

+S (cos - 1) rf t + 2/i /iS' sin cos <p + f *) = (49) 

From Eqs. (14ip and ( 1451) we obtain £* = and r? 2 = rj 2 (z,t,(p). Together 
with Eqs.(T43]) and (gBJ) we get 77 1 = 77 1 (£, 0). From Eq.(g2]) r\\ B = 0, hence 
r/ 1 = a{t)6 + &(£). A substitution of the result into Eq. fj4"TI) yields a = 0, 
b = and r? 2 + £ = 0. This means r? 1 = and = g [Eq.(PJ]. Together 
with Eq. fT44l this produces ry 2 = and = £t = 0, i.e. £* =const and 
£ 2 = £ 2 (^)- Using Eq. fT48l) we get after splitting over variable £ z =const 
and r/ 2 = rf (t). It means that rj 2 , = and we get r? 2 = from Eq. (1491) . By 
collecting all results together we obtain 
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